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Abstract 

In this study we explain the role of applied magnetic field in inelastic conduction properties of a 
Quantum Dot coupled with an oscillator . In the presence of strong applied magnetic field coulomb 
blockade effects become weak due to induced Zeeman splitting in spin degenerate eigen states of 
Quantum Dot. By contacting Quantum Dot by identical metallic leads tunneling rates of spin down 
and spin up electrons between Quantum Dot and electrodes will be symmetric. For symmetric 
tunneling rates of spin down and spin up electrons onto Quantum Dot, first oscillator get excited 
by spin down electrons and then spin up elctrons could excite it further. Where as average energy 
transferred to oscillator coupled with Quantum Dot by spin down electrons will further increase by 
average energy transferred by spin up electrons to oscillator. Here we have also discussed that with 
increasing Quantum Dot and electrodes coupling strength phononic side band peaks start hiding 
up, which happens because with increasing tunneling rates electronic states of Quantum Dot start 
gettting broadened. 
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A. Introduction 



In recent years, much attention has been focused on the concept and realization of nano- 
electromechanical systems (NEMS) [l- 8] as a new generation of quantum electronic devices. 
A large number of new experimental techniques have been developed to fabricate and per- 



form experiments with NEMS in the quantum regime. Examples of high- 



cal nano-structures that have been produced are nano-scale resonators 



quantum dots or single molecules [l ll- If]] . cantilevers 



requency mechani- 



18 



more recently graphene sheets 



20j and carbon nanotubes 



10|, semiconductor 



19 
21 



vibrating crystal b earns |9j, and 



221 ] . These devices are expected 



to open up a number of future applications including nanomechanical transport effects, sig- 
nal processing which could be used in fundamental research and perhaps even form the basis 
for new forms of mechanical computers. Many theoretical methods and models have been 
designed in order to account for the behavior of different types of NEMS system and to 
make predictions and proposals for future experiments. 

In general, there are two different theoretical formulations that can be used to study the 
quantum transport in nanoscopic systems under applied bias. Firstly, a generalized quan- 



tum master eq uation approach 23-32] and secondly, the nonequilibrium Green's function 
formulation 33 3^] . The former leads to a simple rate equation, where the coupling between 
the dot and the electrodes is considered as a weak perturbation and the electron- phonon 
interaction is also considered very weak. In the latter case one can consider weak and in- 
termediate electordes to system and electron-phonon coupling. The nonequilibrium Green's 
function technique is able to deal with a ver y b road variety of physical situations related 



to quantum transport at molecular levels 



36 



37j . It can deal with strong non-equilibrium 



situations and very small to very large applied bias. In the early seventies, the nonequilib- 
rium Green's function approach was applied to mesoscopic transport 38H40| by Caroli et al., 
where they were mainly interested in inelastic transport effects in tunneling through oxide 
barriers. This approach was formulated in an elegant way 4lM43| by Mier et al, where they 
have shown an exact time dependent expression for the non-equilibrium current through 
mesoscopic systems. In this model an interacting and non-interacting mesoscopic system 
was placed between two large semi-infinite leads. In most of the theoretical work on NEMS 
devices since the original proposal, the mechanical degree of freedom has been described 
classically /semiclassically 23|, |44( or quantum mechanically 24j-|26|, |45|, |46| using the quantum 
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master or rate equation approach. In the original proposal, the mechanical part was also 
treated classically, including the damped oscillator, and assuming an incoherent electron tun- 
neling process. This approach is based on a perturbation, weak coupling and large applied 
bias approximations, whereas the Keldysh nonequilibrium Green's function formulation can 
treat the system leads and electron-phonon coupling with strong interactions [47] for both 
small and large applied bias voltage. The transport properties have been described and 
discussed semi-classically/classically but need a complete quantum mechanical description. 
A theory beyond these cases is required in order to further refine experiments to investigate 
quantum transport properties of NEMS devices. In the quantum transport properties of 
these devices; the quantized current can be determined by the frequency of the quantum 
mechanical oscillator, the interplay between the time scales of the electronic and mechanical 
degrees of freedom, and the suppression of stochastic tunneling events due to matching of 
the Fermionic and oscillator properties. 

In the present work, we consider a spin dependent electron transport through a quantum 
dot connected to two identical metallic leads via tunneling junctions. A single nanoelec- 
tromechanical oscillator is coupled with quantum dot and gate voltage is used to tune the 
levels on the dot. The application of strong magnetic field induce Zeeman splitting in spin 
degenrate eigen states of quantum dot. As a result spin down states moves lower and spin 
up states move higher than the degenrate spin eigen states of quantum quantum dot, and 
thus offers the different channels of conductance for spin up and spin down electrons. In 
the presence of strong applied magnetic field the coulomb blockade effects will be weak due 
to Zeeman splitting and we theoratically included it by mean field approximation, which is 
quite resonable approximation for tackling weak interactions. Although electron transport 
through mesoscopic systems in the presence of Zeeman splitting has been an active area of 
research 57H59I] . In our calculation the inclusion of the oscillator is not perturbative which 



enable us to predict strong electron phonons coupling effects in NEMS system. Hence, our 
work provides an exact analytical solution to the current-voltage characteristics, conduc- 
tance, coupling of leads with the system, and it includes both the right and left Fermi- level 

45 0|. where the cou- 
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response regimes. However, we have used wide-band approximation 
pling between leads and quantum dot is taken to be independent of energy. This provides 
a way to perform transient transport calculations from first principles while retaining the 
essential physics of the electronic structure of the quantum dot and the leads. 



3 



B. Model Hamiltonian 

Our mesoscopic system consists of a Quantum Dot(QD) coupled with an Oscillator to 
include the role of phonons effects in conduction through a QD. Application of external 
applied magnetic field induce Zeeman splitting in spin degenerate eigen states of QD. This 
constitute microscopic part of mesoscopic system. To incoporate coulomb blockade effects in 
QD we use mean field approximation which is useful for weak interaction. As in presence 
of strong magnetic field coulomb blockade effects will be small because of Zeeman splitting. 
Hamiltonian of the present microscopic system would be, 



HQD+Osdiiator = ^(eg + -^Bgo .B)d\d a + a a (a + a ] )dld a + u (a f a + -) (1) 

a a 

The first term represents two discrete energy levels in QD, which orginates because of 
magnetic field induce Zeeman splitting, d^, d a (at a) create and annihilate an electron in state 
\a > on the dot (create and annihilate a phonon in state \n > on the oscillator) . Here e a , 
Hb-, 9, B, a a and ui are energy levels of QD electronic-state with spin a ,Bohar magneton, 
Lande g factor, Pauli spin matrix, applied magnetic field, QD and oscillator coupling and 
oscillator vibrations frequency. Second term represents oscillator QD coupling and last term 
represents oscillator energy spectrum. 

In first term of hamiltonian (cr = 1) for spin up electrons and (a — — 1) for spin down 
electrons, 

H = e , e t = e + U (2) 

where U represents coulomb repulsion between spin down and spin up electrons. 

To pass current through this sytem we employ left/right electrodes which constitues 
macroscopic part of of our system. Hamiltonian of left/right electrodes is 

Blends = ^ ] £k,o,v C k,o,v C k,cr,v (3) 

Here represents electrodes electronic states with wave vector k, spin a ,and electrodes 
v (left/right). c\. av {ck )(T)U ) is electron creation (annhilation) operator in electrode v . 

Hopping of electrons between electrodes and QD is defined by the following Hamiltonian, 
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H Hopping = yijTkwdjckw + hermitian - congugation) (4) 

Here T k a v represents electron hopping amplitudes between QD and electrodes. 

We first solve our microscopic system Hamiltonian. Our approach include electron- 
phonon interaction exactly (non-perturbatively). 

To diagonalize microscopic system Hamiltonain , we employ Lang-Firsov 
transformation 51 ] . 



Hqd+Os dilator 

Exp[S}H QD+0s 

dilator 

Exp[S ] ] (5) 

where 

S = J^i^ - a)dld a (6) 
<j 

After diagonalization 



dl = dlExp[— (a f -a) (7) 

d a = d a Exp[ -(a* — a) (8) 

a* = a} -d^dcr ,a = a -S n d a (9) 

Therefore, 



H, 



QD+Osdllator 



^(e ff + — jigger. B - A ff )4<k + w (a f a + -) (10) 



Where ^ = — 

too 



Now the eigen function of the diagonalized Hamiltonian in k-space ( eigen function of 
harmonic oscillator remain same in real and Fourier's space) would be 

1 x - k 2 
\oan >= —=22A n Exp[ — ^-}H n [k a }Exp[-ik a x a }\a > (11) 

v a 

I £2 

\uan >= -f=J2 A n Ex Pl-^Wn[ka]\a > (12) 
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\oan > is state of occupied QD with electron and \uan > is the state of un-occupied QD 
with electron. 

Here x a represents displacement of oscillator due to occupancy of electron in QD. x a = ^ 
and H n [k a ] are usual Hermite polynomials. 

Now the amplitude of the occupied and un-occupied QD electronic state would be, 

Anno- =< oam\uan > (13) 



A mna = 1 == / dk a Exp[-kl]H m [k a ]H n [k a ]Exp[ik a x a ] (14) 

V7r2 m + n+2 n!m! J 



Ann. = V^r^^t-^]^] 1 "-" 11 ^-" 11 ^] (15) 

Here L™ - ™' [^f] represents associated Lagurre's polynomials. 
After diagonalization tunneling Hamiltonian H Hopping will become, 

H Hopping = y~](Tk,a,vdlck,(r, v + hermitian - congregation) (16) 
Where (f h ^ v = T k ^ u Ex P [^ - a)]) 

C. Current from the mesoscopic system 

Current from the (v) electrode to the QD can be calculated by taking time derivative of 
occupation number operator of (v) electrode. 

r(t) = -e(^N v ^ = ie([N v ,H\) (17) 
where (n v = Y, k ,a c l,a,u c k,a,Jj and (h = H QD+0s ciiiator + H Leads + H H opping^ .Therefore, 

r(t) = ie J2 (f k „ (c{ a J„) - Tt„ (#<w)) (18) 



k,0~ 

Now we define electrode and QD coupled lesser Green's function, 
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(20) 



J"(t) = i2elm^f k ^ v G< naM {t,t)\ 



(21) 



To find electrode and QD coupled lesser Green's function we utilize equation of mo- 
tion technique (see 43j, I53M55I] for utilizing equation of motion technique in non-equilibrium 
Green's function theory), 







Q^i e k,a,u ) G unak(J {t t ) T k <JU G un(TUn(T {t t ) 



(22) 



where superscript (t) on Green's function notation represents time ordered Green's func- 
tion. 

Lets define electrode inverse Green's function 



And by using Green's function identity 



.d_ _ 
-% Qt , e k ^ u 



(23) 



(24) 



eq f)22p can be written in the following form, 

" = f Lu I dhGl naiUna (t - ti)dU,(*i - ( 25 ) 
Now by using analytic continuation rule our electrode QD coupled Green's function 



becomes 



56] 



(26) 

We are discussing dc bias situation therefore its useful to work in energy space. Hence 
by using convolution theorem in eq (126]) ,then QD electrodes coupled Green's function will 

be, 
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^Mno-,fccr( e ) ~~ ^k,a,u [^una,uncr{ € ) 9k,cr,u( € ) + ^uncr,uncr( e ) 9k,a,ui € )] (^7) 

Hence current from the mesoscopic system is, 

J" (t) = z2e Im[^ f*,^^ | ^ [G; n<T)WCT (e)^<^(e) + G<„ CT (e)^ )CTiI/ (e)] ] (28) 
Here (g^j,) G r ^ auna represents (electrodes) QD retarded and lesser Green's functions. 



From equation of motion method electrodes lesser and retarded Green's function is given 

by, 

= (e _ e l + —) ; ^( e ) = * 27r5 ( w - ( 29 ) 

and 

Tk,a,vT k(Jll = T). !atl/ T k(TV (30) 

We employ the wide-band approximation, where the energy density of the electrodes is 
taken to be energy independent and — » de ka 



/OO /"OO 
de k «gt*,A e ) = ^ DT ^ T L» / Me - e ka ) de ka = iT v J»{e) (31) 
OO J — OO 

Here (£)) represents constant energy density of electrodes. 

_ <^ W M = OT te „Tt w | _ — ^ = - !f (32) 
By using eqs (1291 - (I32]) in eq (!28|) our current expression, 

r(t) = zeQT / |^ [(Gl„{e) - G a u „{e)) #(e) ~ G< na , w M] \ (33) 

For dc transport current will be uniform J = J L = — J R ,So symmetrize current expres- 

J L - J R 
sion will be J = , 
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J ( £ ) = f E/" ^ K rL ( £ ) - ri? ( e )) G< ( e ) + - G a {e)}{ti{e)T L {e) - f R (e)V R (e)} 

(34) 

Here bold face representations of level width functions T L / R 's and Green's function 
G r ' a,< shows their matrices in microscopic part of the system. 

Now mesoscopic system Green's function could be found from spectral representation of 
Green's function, 

G una( e ) =< una\(e - H QD+ oscillatorI ± iT a )\una > (35) 
r L + T R 

where r CT = - - — — and / = J2oma \ omcr >< oma\ 

nr,a t \ < una\oma >< oma\una > , . 

^ (e - e CT - ^fi B gcrB + A a - (m + ^)u ± iIV) 

El 4 I 2 
/ 1 DiA / i 1 \ _|_ '"pi \ \°' J 

oma ( £ ~ ^ ~ WB9<?B + A CT - (m + 2)^0 ± ^ o" ) 

While doing numerical calculation we have considered oscillator to be at ground state 
so set (n = 0), when spin down elctrons comes onto QD coupled with an oscillator and 
summation (m) plays the role of creations of phonons created by spin down electrons on it. 
While spin up electrons will come on to QD coupled with an excited oscillator at (n = m) 
, excited by spin down electrons. Now summation (m) plays the role of creation of more 
phonons in oscillator by spin up electrons. 



Mesoscopic system lesser Green's function is given by [471] . 



G< na (e) = ^ nCT (e)(/ CT L (e)r^( e ) + f R (e)T R (e))G a una (e) (38) 
Average energy transfer to oscillator by spin up and spin down electrons is defined by, 

E PH = <*> fa"*"* ) (39) 



In eq (l3"9l we have ignored ground state energy contribution to oscillator, which will give 
just a shift in average energy transferred by electrons to oscillator. 

p un „(t,t) = -i2G< na (t,t) (40) 
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0.5 1 1.5 2 2.5 

Applied Bias (meV) 



FIG. 1: Differential conductance plot with respect to applied bias .Here e±(B) = 0.5 meV , 
e t (B) = 1 meV , U = 0.1 meV T = 0.1 meV. 

G< n<7 (t,t) = I ^G< na (e) (41) 

I. RESULTS AND DISCUSSIONS 

A single level QD in the absence of magnetic field has spin degenerate levels. As the 
applied voltage from electrodes become equal to QD energy levels, then spin up and spin 
down electrons start tunneling through the QD. Application of applied magnetic field to 
QD induce Zeeman splitting. As a result spin up level moves higher than unsplitted levels 
and similarly spin down level moves lower than unsplitted levels. Now as applied voltage 
from electrodes is increased first spin down energy level resonates with applied bias and 
then spin up energy level resonates with applied bias. In Fig.(l) we have showed differential 
conductance as a function of applied bias from electrodes. 

Now we explain QD coupling with an oscillator in the presence of magnetic field. At 
zero temperature oscillator will be in ground state. As spin down electron comes onto QD it 
gives energy to oscillator and oscillator moves to excited state. This explain the phononic 
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peaks appearance in differential conductance for spin down electrons. The main peak of 
spin down electron will get shifted due to (A) and amplitude of main peak of spin down 
electron becomes smaller as only spin down electron could give energy to oscillator in zero 
temperature. We have assumed strong dissipation effects with enviornment, which means as 
spin up or spin down electron leaves oscillator it comes to ground state. This rules out accu- 
mulation of energy in oscillator. When applied bias resonates with (e| + (m + |) + U — A) 
then spin up electron channel too get activated along with spin down electron channel. Here 
(m)represents number of phonons produced by spin down electron, and moreover spin up 
and spin down electrons have symmetric coupling with electrodes and QD (Dj- = V J which 
means both spin up and spin down electrons comes onto QD and leaves the QD in the same 
time, and this is quite resonable as for identical metallic electrodes QD electrodes coupling 



strength will be same for 
using ferromagnetic leads 



Doth spin up and spin down electrons, This could be changed by 



521 ] . Therefore spin up electron excites oscillator to even more 



excited state. So we get satellite peaks in differential conductance of spin up electrons. Here 
phononic peaks of spin up and spin down electrons is not same. This happens as excited 
states of occupied QD coupled with oscillator eigen states are not same for different values 
of excitation. See Fig(2). 

Effects of electrodes-QD coupling (Ts) in the presence of applied magnetic field and 
oscillator-QD coupling is of particular importance. As the tunneling rates of spin up and 
spin down electrons from electrodes to QD is increased then energy states of the QD gets 
broadened. In Fig. (3) we have showed that for a fixed value of applied magnetic field and 
QD-oscillator coupling when electrodes-QD coupling are small than Zeeman splitted peaks 
and phononic side band peaks are clearly visible. But as we increase electrodes-QD coupling 
then phononic side band peaks starts disappearing. 

Average energy transferred to the oscillator by spin down and spin up electrons is shown 
in fig(4). Here we could see that spin down electron curve lies lower than spin up curve, 
where as small steps are signature of phonons creation in averge energy versus applied bias 
plot. Spin up electron starts contributing to increase average energy of oscillator where spin 
down electrons ends up, and applied bias resonates with (e^ + (m + |) + U — A J . 
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FIG. 2: Differential conductance plot with respect to applied bias .Here ej.(-B) = 0.5 meV , 
e t (B) = 1 meV , U = 0.1 meF ,r = 0.3 meV,a = 0.37 meF ,w = 0.29 meV. 
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FIG. 3: Differential conductance plot with respect to applied bias .Here ej.(-B) = 0.5 meV , 
e t (B) = 1 meV , U = 0.1 meF ,q = 0.37 meV ,w = 0.29 meV. 
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FIG. 4: Average energy plot with respect to applied bias .Here e±{B) = 0.5 meV , e^-(B) = 1 meV 
,U = 0.1 meV ,a = 0.45 meV ,w = 0.35 meV. 



A. Conclusion 

In this work we have studied inelastic electron transport through QD coupled with an 
nanomechancial oscillator in the presence of strong applied magnetic field. We have ex- 
plained first spin down elctron creates phonon in QD and then spin up start creating phonons 
on it. Due to creation of phonons small steps are produced in average energy transferred to 
oscillator. With increasing electrodes QD coupling strength phononic satellite peaks starts 
hiding up. 
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